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I. INTRODUCTION 



I cannot claim the privilege of having known Professor John Bell as intimately or for as 
long as some of the other participants in this conference to honour his memory. But there 
was a period of about 3 years from 1982 to 1985 when we interacted a fair amount and 
conducted (mostly at long distance) a collaboration trying to understand the then recently 
discovered phenomenon of fractional charge. During this period we had the opportunity of 
hosting him at my Institute in Bangalore. He came home for dinner and met my family. I 
in turn went to CERN , first to write up our work and later for a long visit, when I had 
the pleasure of meeting his wife Dr. Mary Bell and was the beneficiary of their kindness in 
many ways. 

During this interaction I developed , like so many people before me, great respect for 
Professor Bell, not just for his great originality and prowess as a theoretical physicist, but 
equally for the precision of thought and language that he demanded of himself and of those 
who were fortunate enough to work with him. John Bell was more than an intellectual force 
in physics — his was a moral presence. One can ill-afford the loss of such a moral presence 
in today's world, and it is a great tragedy that he was taken away from us when we continue 
to need him so much. 

Long before I met John, I had of course heard of him and his many important contribu- 
tions. In fact the very first time I had to study his work was in the early 'sixties, when I was 
a foot-soldier in Professor Hans Bethe's army attacking the nuclear matter problem. All of 
you know that John had worked on a wide range of subjects, but not all might know that one 
of his early interests was nuclear matter, on which he had co-authored a long review article 
with E.J. Squires way back in 1961. A few years after that as I gravitated towards particle 
physics and quantum field theory I had to study in great detail his landmark discovery with 
Roman Jackiw of the axial current anomaly. 

Despite such overlap of interests, I met him only much later in life when he spent a 
couple of weeks at the Indian Institute of Science in Bangalore where I was then working. A 
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colleague introduced us and he made the usual polite enquiry - "What are you working on 
these days ?" I thought that rather than inflict on him the different problems I was working 
on at that time, I would instead get his "take" on the very interesting but puzzling claims 
of fractionally charged states that had recently emerged. 

As it turned out, he had not heard about these developments. "You mean fractionally 
charged like quarks ? Isn't that by now old and well established ? ", says John , a slight 
note of wariness creeping into his voice. " No, no ", I hasten to add, "I don't mean quarks . 
By fractionally charged I refer to states with fractional eigenvalues of the Number Operator. 
These have been discovered not just theoretically in soliton states of field theory models, 
but are claimed to be present in down to earth systems like polymer chains . What I don't 
understand is how a real polymer made of some finite number of electrons could, no matter 
how it twists itself into a soliton configuration, carry a fractional number of electrons." This 
produced a flicker of a smile from him and I knew I had caught his interest. "Would you like 
to tell me about it ?" , he asked. So I ushered him into my office and closed the door . This 
was a lovely opportunity , having trapped the distinguished John Bell in one's office with 
all avenues of escape sealed off, to summarize for him ab initio the curious phenomenon of 
fractional charge and get his opinion. 

Let me begin the next section with a pedagogical expansion of what I described to John 
in the fading light of that winter afternoon in Bangalore. 

II. FRACTIONAL CHARGE IN FIELD THEORY 

The phenomenon of fractional charge was first discovered by Jackiw and Rebbi in their 
pioneering work in 1976 in the context of soliton states in quantum field theory |J. Three 
years after the Jackiw-Rebbi work, but quite independently and in an entirely different 
context, Su, Schrieffer and Heeger [[|, argued that a similar phenomenon can occur for 
soliton states in the long chain molecule trans- Polyacetylene. See also the work by Rice j3]. 
Subsequently Jackiw and Schrieffer [Q wrote a paper drawing attention to the similarity 
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between these results discovered in entirely different areas of physics. 

Let us begin with the Jackiw-Rebbi work . They discussed the phenomenon in both 
one- dimensional and three-dimensional models. Let me for simplicity discuss just the one- 
dimensional example, which contains all the essential ideas. Consider in (1+1) dimensions 
a Fermi field ty(x,t) coupled to a scalar field <£>(x, t) through a Lagrangian density 
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and 



C F = ^[id^ - m$(x,f)J# (3) 

This example corresponds to a quartic double well potential in $. It is straighforward to gen- 
eralise the results to other potentials U(<&) with symmetric double wells. The phenomenon 
of fractional charge occurs in the soliton sector of this system. To a wider readership it may 
be helpful to recall what is meant by the soliton sector and how fermi fields are treated in 
such sectors ( [f|). 

The bosonic sub-system Cb, in the absence of the Fermi field, has as its field equation 

(|r - v 2 )^ - + 3 = o (4) 

This has the following two lowest energy classical solutions : 

<p(x,t) = ±l (5) 

As is well known the presence of these two degenerate classical solutions indicates sponta- 
neous breaking of the $ <-> — $ symmetry of the Lagrangian. Around each of these classical 
solutions a separate vacuum state and a whole tower of Fock states can be built correspond- 
ing to the two phases this system permits. We will call these the vacuum sectors. 
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In addition this system also has two other static (time independent) solutions. Those 
are the topological soliton solutions - the so called kink solution and its reflected partner, 
the anti-kink, given respectively by 



As per the general theory of semi classical quantization of quantum fields, one can build two 
other separate towers of states, one around each of these soliton solutions. The topological 
index n = 0(oo) — </>(— oo) becomes a superselection quantum number upon quantization 
and prevents any states from these soliton sectors from decaying into the vacuum sector. In 
short, we have four sectors of states for this bosonic system : two are the vacuum sectors 
built around <fi = ±1 and the other two are the soliton and antisoliton sectors built around 
the kink and antikink solutions in eq (|G|). For more details on soliton sectors of states and 
their properties see ( 0). 

Now consider the full system in ([]]) including the Fermi field. 



To leading order in % one can replace in each sector the Bose field operator <3> occurring 
in the fermi field lagrangian fl3|) by the corresponding classical solution . Thus in the vacuum 
sector built around 0=1, the fermi lagrangian @) reduces in leading order to 



This is just the free Dirac system with mass m discussed in textbooks. In order to contrast 
with what later happens in the soliton sector , let us recall why the Number (Total charge) 
operator has only integral eigenvalues for the free Dirac system. Let us denote by Uk{x) and 
Uk(x) the positive and negative energy spinorial solutions of the Dirac equation 




(6) 



A. Vacuum sector 




(7) 



(— iad x + (3m)uk{x) 



E k u k (x) 




(8) 
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where E k = +\/k 2 + m 2 and spinor indices have been suppressed. The Dirac matrices for 
this 1+1 dimensional system can be taken to be a = o<i and (3 = <j\. The Dirac field is 
expanded in terms of these solutions and the destruction operators b k and d k obeying the 
usual anticommutation rules . 

*(x,t) = E [b k u k e- iEkt + d\u k e iE ^] (9) 

k 

The vacuum state in the 0=1 sector is given by the conditions 

bk\vac) = dk\vac) = (10) 

with all the bosonic oscillators being in the ground state. 

Note that the third Pauli matrix 03 acts as the charge conjugation matrix. It anti 
commutes with the Dirac hamiltonian in and generates for every positive energy solution 
Uk(x) of energy E k the corresponding negative energy solution of energy — E k : 

o?,u k {x) = u k {x) (11) 

Hence all modes of the expansion ( |9|) come in pairs with positive and negative energy. There 
are no zero energy solutions in the free massive Dirac equation. 

Finally consider the charge density operator (which is really the number density operator) 

1 



&{x,t) ,^{x,t) 



(12) 



This commutator form of p is standard; it is designed to be regularised and odd under charge 
conjugation . Inserting the mode expansion @ and using the orthonormality of the Dirac 
solutions the total charge becomes 

Q = J dx p(x, t) 

= \ E A] + [4,4]) 

= £ ( {b\ b k - 1/2) - (44 - 1/2)" 

= E U h - d\d k ) (13) 



Notice that the half-integers cancel term by term because of the existence of paired positive 
and energy modes. Since 

(&I h) 2 = b\ b k , (4 d k f = 4 d k (i4) 

these have eigenvalues of only or 1. Hence the familiar result in the vacuum sector that 
the charge operator has only integer eigenvalues. 

B. soliton sector 

Let us repeat exactly the same steps in the soliton sector. Now we have to substitute for 
the Bose field operator in the Dirac Lagrangian the classical kink function in eq (|j). The 
corresponding Dirac equation now becomes 

x 

(—iad x + (3 m tanh—j= )ip(x) = Eip{x) (15) 

For the two components ^1,2 of the Dirac spinor, this yields the coupled equations 

x 

( —d x + m tanh—i= )ip2 — E 
v 2 

x 

( d x + m tanh—i= = E ip 2 (16) 
v2 

This equation will also have a set of positive energy solutions r] k (x) with some associated 
energy E k . We need not find these solutions explicitly. But since the charge conjugation 
matrix 03 again anticommutes with the Dirac hamiltonian, we know that for every positive 
energy solution r)k{x) there will exist a negative energy solution fjk(x) with energy —E k . But 
now there is also an unpaired zero-energy solution 



(A expi — m J x dy tank (y/\/2) 
' 1 <17) 

Such a normalisable solution to (16) exists because the soliton function tanh^ which forms 
the background potential for the Dirac spinor tends to opposing limits ±1 as x — >• ±00. In 
infinite spatial volume this solution has no partner. It is self charge conjugate : 03770 = t]q. 



The mode expansion of the Fermi field operator now becomes 

tt(z,t) = E [hVk(x)e' iEkt + 4%(x)e i£fc '] + ar, (x) (18) 

where a is the destruction operator for the zero-mode. 

Unlike the vacuum sector built around 0=1 which had a unique ground state ("the 
vacuum" in that sector), in the soliton sector there are two degenerate ground states because 
of the existence of the fermionic zero-mode. They are \sol) and \sol) obeying 

a\sol) = bj~\sol) = dk\sol) = (19) 

and 

\sol) = at] sol) ; a \sol) = \sol) (20) 

These are the two basic quantum soliton states of this system. They are energetically 
degenerate, but are distinguishable by their Charge. 
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Notice that the piece (-1/2) coming from the zero- mode commutator remains uncancelled 
because it does not have a charge conjugate partner. Since the operators b\ b k , d\ c4 and 
a) a all have eigenvalues of or 1, it then follows that the total Charge (Number) operator Q 



has half-integral eigenvalues . It should be emphasized that eq([21|) is an operator equation 



for Q. The half-integer appearing in it will be reflected in its eigenvalues and not just its 
expectation values. In particular, the two degenerate soliton states have eigenvalues of ±1/2 
respectively for the total number operator Q : 

Q \sol) = -(l/2)\sol) Q \sol) = (1/2) \sol) (22) 



This, in its barest form, is the prototype example of the Jackiw-Rebbi discovery that 
in some field theories there can be states carrying fractional eigenvalues of the Number 
operator. For a fuller discussion including similar results in three-dimensional models see 
their original paper [I]. 

III. POLYACETYLENE 

This discovery by Jackiw and Rebbi was clearly very remarkable. But it is not easy 
to "understand" the result physically in terms of our familiar intuition with quantum field 
theory. How can the Number Operator which is widely used to count the number of particles 
in QED and so many other theories, and which even in the Higgs model used by Jackiw and 
Rebbi has only integer eigenvalues in the vacuum sector, yield fractional values in another 
sector ? Yet the proof given is so simple and transparent that, stare at it as we may, we 
have no choice but to accept the result within the parameters of its derivation. In my own 
attempts to make peace with the result (before embarking on the more detailed study with 
John), I loosely attributed it to the vagaries of the infinite degrees of freedom of continuum 
field theory. [ We know from high school that formally summing divergent series of integers 
can yield fractions.] But clearly some further clarification of this phenomenon was called 
for. 

The need for clarification became more compelling when a similar result was derived not 
in a model field theory, but in a down-to-earth experimentally accessible polymer system by 
Su Schrieffer and Heeger ||. Three years after the Jackiw-Rebbi work, but quite indepen- 
dently they showed that the same phenomenon of fermion number 1/2 can occur for soliton 
states in the long chain molecule trens-Polyacetylene. 

That the same phenomenon occurs in Polyacetylene (we will henceforth drop the prefix 
"trans") as in the Jackiw-Rebbi field theory is not accidental. As pointed out by Jackiw and 
Schrieffer , the former system has the same structure in the continuum limit as the field 
theory model. Polyacetylene is the molecule (CH) n with large n, where the Carbon ions 
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form a long chain with the H atoms sticking out transverse to the chain. If we consider the 
flow of electrons along the chain, the system can be viewed as consisting of electrons and 
bosons (the phonons of lattice vibration of the Carbon ions along the chain) in one space 
dimension just as in the field theory model . 

The Hamiltonian for the electron-lattice system can be taken, for each of the two spin 
states of the electron and for small ion displacements as (see [§J]) : 



H = E 



2 

+ K/2 K +1 - u n ) 2 ) + (D ] n+1 D n + h.c.)(u ll + l - u„ 1/(2*) 



(23) 



where n labels lattice sites, p n and u n are ion momenta and displacements respectively , D n 
is electron destruction operator and K and a are constants. All other constants have been 
absorbed into definitions for simplicity. 

Our primary interest is in the fermions . As far as the lattice vibrations go, let us just 
state the Su et al result that the lattice system undergoes dimerisation doubling its 
spatial period (a Peierls transition) . There are two degenerate ground states. In one , the 
mean value of the displacement u n instead of being zero takes the staggered value of 

K = 4(-l)"u„ = 1 (24) 

for all n. In the other 

<j> n = -1 (25) 

The two degenerate ground states further provide the possibility of domain wall (soliton) 
configurations connecting the two phases , i.e. <p n — > ±1 (or =p 1), as n — > ±oo respectively. 
Thus the boson subsystem has four sectors of states - two of them being soliton sectors - 
just as in the field theory discussed in the previous section. 

Furthermore, the boson coordinates u n again act as the background potential for the 
electrons . In fact the electronic part of the hamiltonian ( |2"3"D becomes in the continuum 
approximation exactly the same as that of the Dirac field theory. John and I , as part of 
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our work on these systems ||, offered a simple derivation showing that the polyacetylene 
hamiltonian is just a realisation of the Kogut-Susskind lattice regularisation of the Dirac 
system H. Let me sketch that derivation. The electronic part of (E3f) is 



H. 



dec 



E 



Di +l D n (u n+ i - u n - 1/ (2a) 



+ herm.conj 



(26) 



Define staggered variables 



B 2r -1 = ( — I) 7 " D 2r -l ', C 2r = ( — 1)' D 2r 



with <f) n = 4(— l) n u n for all n as already defined above. Then, 



(27) 



H, 



dec 



E 



D 2r+1 D 2r (u 2r+ i - u 2r - l/2a) + D\ r D 2r _ x (u 2r - u 2r -i - l/ 2a )j + h.c 



(1/4) Bl +1 C 2r (0 2r+1 + 2r + 2/a) + C^5 2r _! 



^2r — <+>2r~l 



2/a) 



r L J 

= E( ^ r+1 2fl ^ r " 1 ) c '2r + (i/4)i4 +1 c 2r . (0 2r+1 + 2r ) + (i/4) J B| r ._ 1 a 2r (0 2r + 2r _x) 

+ /i.e. (28) 

In the continuum limit, as a — > and 2 r ~~ * J ^ this reduces to 

(dBt „ „ + <9 B 



H, 



dec 



dx 



Kax C + C* w ) + M&C + CtB) 
j dx (- Bt*Z + Ct^) + 0(x) (StC + CtS ) 



(29) 



This is just the Dirac Hamiltonian with Yukawa coupling : 







(30) 



with the upper and lower components of the 2-spinor being identified respectively with the 
odd-and even- site electron operators : 

'B(x) 
C(x) 

The charge operator operator at each site also has the same commutator form as in the 
continuum theory's charge density : 



(31) 
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Pn = D\D n - (1/2) = (1/2) 



Dl,D n 



(32) 



where the 1/2 subtracted at each site can be attributed to the neutralizing background 
(ionic) charge per spin state. 

Given this mapping from the polyacetylene system to the field theory model used by 
Jackiw and Rebbi, it is not surprising that fractional charged states arise in the former 
system too. 

But in the context of a real polymer the result becomes still more mysterious. A molecule 
of Polyacetylene has after all some finite number of electrons and it is hard to imagine how 
, even if it twists itself into some topological soliton state, it could have a fractional number 
of electrons. The need to understand this better becameall the more compelling because 
Su, Schrieffer- and Heeger referred to experimental signals supporting the result . See for 
instance ref ( [0]). [ Because of the presence of the two spin degrees of freedom in actual 
polyacetylene, as distinct from a truly 1+1 dimensional model, the fractional charge gets 
doubled and becomes integral; but a signature of the effect can still be observed in the form 
of an unfamiliar charge-spin combination of excitations. There will be neutral excitations 
with spin and charged excitation which are spinless.] 

It should be emphasised that an expectation value of one-half would cause no surprise 
. Obviously that can easily arise for an operator with integral eigenvalues. For instance if 
we take any particle in the ground state of a double-well potential and look for it in one of 
the wells, you will find it there half the times and not find it there the other half, giving an 
expectation value of 1/2. While in the Su et al paper, fractional values were not claimed as 
eigenvalues, in the prototype field theoretic case we have explicitly seen in sec. II that the 
eigenvalues themselves are half-integers. This is what puzzled me and I was fortunate in 
having John Bell join me in worrying about it further. 
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IV. EIGENVALUE OR EXPECTATION VALUE ? 



Now let me summarize what John and I did in our attempt to understand better what 
these half integral values of the number operator mean [f§, . [ We later learnt that 
Kivelson and Schrieffer had also been independently investigating the eigenvalue status of 
fractional charge around the same time with essentially a similar resolution of the issues |TU 



See also the work of Jackiw etal which provided further clarification [II]. I will describe here 
only the work done with John Bell.] Our strategy was to treat the infrared and ultraviolet 
limits of the problem more carefully. In one paper ||, we looked at the continuum field 
theory used by Jackiw and Rebbi, but in a finite volume, which could be later taken to 
infinity. In a subsequent paper we considered what would happen if we put an ultraviolet 
cut-off as well, by going to the finite polyacetylene chain |§. 

Let me begin with the field theory example of Jackiw and Rebbi which was discussed in 
section II. The result derived there that the total charge operator has half-integral eigenvalues 
cannot be disputed in the infinite spatial volume case (L — > oo in the one dimensional field 
theory model considered there). However suppose we kept the spatial volume 2L finite to 
start with, and later go the L — > oo limit ||. In particular , in the soliton sector, we go 
back and solve the same Dirac equation 

x 

( —d x + mtanh^= )ip 2 — E ipi 
x 

( d x + m tanh—i= = E %p 2 (33) 
v2 

but in a finite volume (— L < x < L) with appropriate boundary conditions on the spinor tp a 
at x = ±L. The natural boundary conditions which leave the Dirac hamiltonian hermitian 
are ipi^{—L) = ?/> lj2 (L). Once again every positive energy solution 77^ will have a negative 
energy partner fjk- We do not need their details. But now there are two zero energy solutions 



'A expi — m Jn dy tanh (y/\/2) 
//«.(• r) = ( V q 1 1 (31) 
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and 



/>.(.<■) = ( / \ ) (35) 

A ea;p( — raL + m J dy tanh (y/y/2) 



While the solution (|34|) is the same one as before , localised around the origin, the second 
solution (|35| ) has support near the edges ±L. In section II where we started with an infinite 
volume problem, this second solution could not be entertained. But for any finite L however 
large this second solution, finite everywhere and normalisable, is certainly present. Each of 
these two zero-modes is charge self-conjugate and both will appear in the expansion of the 
Dirac field operator, along with the non-zero energy modes : 

*(x,t) = £ [bkVk{x)e- iEkt + d{fj k {x)e iEkt } + arj (x) + c ] fj (x) (36) 

Correspondingly there will now be four degenerate ground states in the soliton sector, 
namely, \sol) defined by 

bk\sol) = dk\sol) = a\sol) = c\sol) = (37) 



and \sol) = a^sol) , \sol) = c^sol), and \sol') = a^c^\sol). 

The expansion ( p6] ) when inserted into the definition of the charge operator 

1 r L 
Q = - dx 

V 2 J-L 



(38) 



will yield 

Q = E (( b kh - 1/2) - (44- 1/2)) + (ot - 1/2) + (etc - 1/2) (39) 



k x 

We see that there are now no unpaired terms of 1/2 and the total charge operator has 
only integral eigenvalues. One can now let the volume 2L become larger and larger and 
at each stage both zero modes will be present and the total charge will have only integral 
eigenvalues. 

That one can restore integral valued charges by this strategy of starting with a finite 
volume and then letting 2L — > oo can also be understood in a different way. Note that if 
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we had closed our line (— L < x > L) by identifying the points x = ±L, that is equivalent 
to putting the system on a circle with a soliton at x=0 and an antisoliton at x = L. This in 
turn amounts to working in the vacuum sector , where we saw that the charge is integral. 

While integral charge has been thus restored, if we were to stop here we would have lost 
all the important physics unearthed by Jackiw and Rebbi. Notice that while there are now 
two zero modes, both localised with width of order 1/m, one of them, fj Q (x) given in eq.([35"|) 
is stuck at the edges of the sample. But the other zero mode r]o(x) in (HD > is m the middle 
and can move about as an excitation, carrying charge 1/2 in some sense that needs to be 
made more precise. 

To pick just the charge carried by this middle zero-mode, we can try and define a "partial 
charge operator" : 



where l,L — > oo in the order L » I » 1/m. In this limit the partial charge Pi clearly 
includes the central zero mode while excluding the second zero mode at the edges ±L. 
However, the soliton ground state will not be an eigenstate of Pi, but only of the total 
charge operator Q in eq(|39|). The charge density operator p(x) given in eq(|12"D and partial 
charge operators such as Pi will excite particle-hole pairs when acting on the soliton state 
| sol). But consider expectation values in this soliton state. 



where in the last line, contributions from non-zero modes cancel for each k by charge con- 
jugation. Hence 




(40) 



( sol\ p(x) \sol) = (-1/2)53 Vk( x ) Vk{x) - fj k * {x)fj k (x) 




(41) 




(42) 
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Because this is not an eigenstate of Pi there will be fluctuations. Since 

(P, - <Pi >)\sol) = ]T f dx ( Vk (x) f, k ,{x)) \k,k') (43) 

k,k' 1 

We have 

(sol\ (Pi - <Pi >) 2 \sol) = £| fdx r)k(x) f,k>{x) I 2 (44) 

k,k' 1 

The integrand is positive and hence the fluctuations are non-zero. In fact they diverge 
logarithmically. 

Such large fluctuations in partial charges are not special to soliton states. They occur 
even in familiar vacuum states and are because of the sharp boundary of the defining do- 
main. One can kill the fluctuations while still retaining a value of 1/2 by defining a more 
sophisticated partial charge operator, one with fuzzy edges. Let 

I rl+d 

Ifc-jjf dtP„ (45) 

When L,l,d — > oo with L >> l,d then once again, 

( sol\ P l4 \sol) = -(1/2) (46) 

To obtain the fluctuations of this operator note that [p(x) — < p(x) >} connects the state 
| sol) to the particle- hole states \k, k') . We therefore have 

(k,k'\Pi - <Pi>\sol) = (k,k'\J dx (p - <p>) \sol) 

— ~r=, 7r~:(k,k'\ [ dx—r\sol) 

%{E k + E k ,y ' 1 J -i dV ' 

' (k,k'\ - \sol) (47) 



i(E k + E k 

by current conservation. Note that the current j has matrix elements (k,k'\j(x)\sol) = 
r)l(x)iafj k '(x). Hence, 

1 r l+d 1 



(48) 
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To obtain the exact result for this one must know analytic expressions for all the non-zero 
energy modes rjk{x) and fjk(x). We had in fact obtained these solutions for the simplified 
case where the soliton function is taken as a step function <p{x) = Q(x) instead of tank x. 
This simplification does not affect the issues of interest to us. But even quite generally one 
can see that for large I far away from the soliton center (I » 1/m) the solutions %(/) will 
be trigonometric functions of kl. Since Ef. ~ k for large k , the matrix element in ( [47]) will 
behave as 1/k for large k and give rise to the logaritmically divergent fluctuations in eq. 
(fy|). But the charge operator with fuzzy edges P\^ defined in ([45]) has an additional integral 
over dl' which will bring down an additional convergence factor of 1/k . There is also the 
extra factor 1/d in the definition which goes to zero. Altogether 



where the function f(k, k')) has a convergent sum over k, k' . Clearly this vanishes as d — > oo. 
In other words, one can define a suitable partial charge operator with fuzzy edges for which 
the soliton state has a value of -1/2 with no fluctuations ! This fractional charge then is an 
eigenvalue. 

It is very possible that if one were to measure the charge in the neighborhood of the 
soliton's center, one is using some such partial charge operator with fuzzy edges. In that 
case the measured value of 1/2 can be promoted to the status of an eigenvalue. Further, 
since these fractions are eigenvalues of complicated partial charge operators and not of the 
total charge (Number) operator Q, there is no conceptual problem reconciling the result 
with our intuition about the latter. 



{sol\(P l4 



<P lid >) 2 \sol) = J2 \(k,k'\P ltd - <Pi, d > \sol)\ 2 




(49) 



A. Polyacetylene re-visited 



Let us describe briefly what a corresponding analysis |§ yields for polyacetylene which 
is a lattice chain of some finite number N of sites amongst which the electrons hop. I will 
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emphasize just the additional insights that emerge here due to the fact that it is an even 
better regulated system than the field theory in finite volume that we just studied. It has 
only a finite number of degrees of freedom and hence no ultraviolet or infrared problems. The 
lattice Dirac Hamiltonian in eqfl25|) yields the following coupled equations for the electron 
wave functions b n and c n at odd and even sites respectively : 

E b n = (l/4)c n+1 (0 n+1 + n - 2/a) + (l/4)c n _i (0 n + n _ x + 2/a) 

E c n = (l/4)6 n+ i (0 n+ i + 0„ + 2/a) + (l/4)6 n _! (0 n + (j) n _ x - 2/a) (50) 

Clearly, given a solution for some positive energy E we can get another solution of negative 
energy (-E) by replacing b n — > —b n , c n — > c n , E —> —E in the above equations. Thus 
positive and negative energy solutions again come in pairs. But suppose the lattice had 
altogether only an odd number N of sites. Then the Dirac operator can have only the same 
odd number of independent solutions. Therefore there must be one or more zero-energy 
solutions. 

Notice that this argument for the existence of zero modes has nothing to do with whether 
the background boson field is in the topological soliton sector or the normal vacuum sector. 
All you need is an odd number of lattice sites ! One can explicitly verify this result. Take 
N = 2M + 1 with, say, M odd. A zero energy solution of the lattice Dirac equation ( |50| ) 
clearly exists in which the even site variable c n = for all n while the odd site variables b n 
satisfy 

b n +i 2 - q(0 n + n _i) 
K-i 2 + a{<f) n + (f) n+ x) 

This will be true for any background Bose field <p n . In the uniform phase (the "vacuum 
sector here) <p n = 1 for all n. Hence the solution is 

b n = fro exp(-Kan) ; c n = (52) 

where the constant K obeys tanh(Ka) = a. If <ft n is a soliton , say having the simple form 
<fr n = A then the fermion zero mode is 
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b n = b exp(— Ka\n\) ; c n = (53) 

Given the single zero mode in both sectors, one can show that the ground states in both 
sectors will have total charge of (-1/2) as eigenvalues. By simply adding + 1/2 to the 
definition of charge one can trivially make eigenvalues integral in both sectors as was the 
case in the finite-volume field theory. Similarly, if the number of lattice points N had been 
even, there would have been either no zero modes or an even number of them in both sectors, 
leading to integral total charge. 

Thus the important lesson we learn is that as far as the total charge is concerned, there 
is no half integral eigenvalue in the soliton sector as compared to the uniform phase. This is 
true in finite chains even if there is a single isolated zero mode in the soliton sector, because 
a zero mode will then exist in the vacuum sector too. 

Of course one can also define for lattice chains partial charges that pick up the zero 
mode in the soliton sector but not the one in the vacuum sector. [ Recall in the example 
mentioned that the zero mode in the soliton sector given in (|52]) is located near the center 
whereas in the vacuum sector even though a zero mode ([53]) was present it was near one of 
the edges.] As happened in the field theory case discussed earlier in this section, such partial 
charge operators would have a half-integral value in the soliton state. This value would get 
promoted to the status of an eigenvalue with no fluctuations if, in the limit of a very long 
chain, the operator covered a region with suitably fuzzy edges (see ref || for more details). 

V. CONCLUSION 

I have tried to outline in the last section the work that John Bell and I had done on 
Fractional charge, preceded in earlier sections by some background on where matters stood 
when we began our work. I believe our work had thrown light on which operators do have 
eigenvalues of 1/2 (as distinct from expectation values ) and which don't. In the process we 
had also laid to rest the concern that had motivated us to study this problem — which was 
that the total number operator of a physical system should not have fractional eigenvalues. 
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That would have made no sense. We showed that the total charge operator defined to 
have integral eigenvalues in the absence of solitons will, in a well regulated theory, continue 
to have integer eigenvalues even in the presence of solitons. But suitable localised partial 
charge operators can have fractional eigenvalues and fractional charge should be interpreted 
as corresponding such an operator. Indeed these may well be what charge measurements 
in experiments employ. As mentioned already, Kivelson and Schrieffer [|l(| had come to the 
same conclusion independently. 

The work discussed above is nearly 20 years old. This article is primarily devoted to 
John Bell and his work , so this is not the place to discuss in any detail other developments 
on this topic around that time or since then. However, for the sake of completeness let me 
summarize in a few sentences some of these developments . 

Soon after charge 1/2 states were discovered in field theory and polyacetylene, the pos- 
sibilities of charge at other fractions were unearthed both in polymer physics and in model 
field theories ( see for example [12j and |13[) . 



A major addition to the list of fractionally charged objects was the theoretical discovery 
by Laughlin in an entirely new arena. He showed that quasi-particles in fractional quan- 
tum Hall (FQH) systems carry fractional charge [TJJ]. Subsequently these have also been 



experimentally observed in "shot-noise" experiments To the best of my knowledge 

the quantum mechanical status of this quasi-particle charge in FQH has not been ana- 
lyzed as carefully as in the case of polyacetylene or field theory models. The former have 
been identified through ingenious but indirect arguments involving plasma analogies and 
Aharanov-Bohm/ Berry phases rather than by an explicit analysis of the charge operators 
and their eigenstates. The FQH quasi-particles are more difficult states to study. They 
correspond to genuinely correlated many-body states as distinct from the polyacetylene and 
Jackiw-Rebbi soliton states considered above, which could be analysed in terms of single 
particle solutions of the Dirac equation. But broadly speaking, the status of the fractional 
charge in FQH is believed to be similar to those discussed above. It is localized near the 
center of the quasi-particle wavefunction and the complementary missing fraction is believed 
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to be near the edge of the Hall sample. 



Finally, very recently, it has been claimed by Maris ||T6[ that some "bubbles" formed 
in liquid helium which contain a trapped electron can then fission into daughter bubbles 
each of which carries a fragment of the original electron . However, in an analysis of this 
phenomenon, the pioneers of fractional charge Jackiw, Rebbi and Schrieffer have argued 
that these fractions in helium bubbles are just expectation values of the sort familiar in 



local charge measurements in any double well quantum system ||17|| . This is an ongoing area 
of work. 
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